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SYNOPSIS 


The problem of stability in a finitely deformed 
solid can be studied either by using the method of 
superposing small peturbations on large deformations or, 
by employing the stability criterion with the help of 
variational technique. The thesis presents a detailed 
study of the following four cases: 

(i) Instabilities in a rectangular viscoelastic 
solid under axial load in a state of plane 
strain, 

(ii) Axisymmetric instabilities in a viscoelastic 
circular cylinder under axial loading. 

(iii) Stability of a thin-walled viscoelastic 
cylinder under finite twist. 

(iv) Stability of a thin rectangular viscoelastic 
plate under in-plane loading. 

The chapter-wise break-up of the thesis is as 

follows; 

Chapter one presents a brief literature survey in 
the subject of viscoelasticity and stability in a visco- 
elastic continuum. It includes a complete formulation of 
the rate boundary value problem and presents the derivation 
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of the stability criterion based on the Kelvin-Dritchlet 
concept of stability of a system in motion. Starting 
from a general viscoelastic material, the relevant 
constitutive relations are obtained in the rate form and 
are finally specialized for the case of an incompressible 
viscoelastic material which this thesis is chiefly 
concerned wibh. 

In chapter two, the sbability phenomenon of a 
rectangular viscoelastic solid under plane strain in 
compress ion/ tension is analysed. By solving the governing 
differential equation of motion and studying the variation 
of the resulting displacement wibh time, the critical 
conditions are obtained for instability, A detailed indepbh 
sbudy of the behaviour of the resulting characteristic 
equabions is presented which makes the discussion on the 
stability phenomenon more rabional and complete. It is 
found that the instability in the symmetric mode is ruled 
out complobely. The results for a corresponding elastic 
material arc obtained as a special case ; these results 
are similar to those obtained by earlier researchers but the 
manner in which they are obtained is mathematically (and 
also logically) more rigorous. 



stability of a viscoelastic solid, cylinder of 
f ini be dimensions under axial load is examined in chapter 
three. Using the formulation developed in chapter one 
and following the approach adopted by Cheng et al, ('1971) 
and Kumar and Niyogi (1982) for an elastic material, a 
condition is derived for the existence of a non-trivial 
solution. Only the axisymmetric mode of instability is 
considered. The instability condition is numerically 
solved to obtain the critical stress for different values 
of cylinder dimensions. 

Chapter four is devoted to the investigation of 
stability of a thin-walled long circular viscoelastic 
cylindrical shell under finite twist. The stability 
criterion presented in chapter one is used to derive the 
condition for stability. The approach is similar to that 
of Neale (1973) and is enlarged to include the case of 
viscoelastic materials. By studying the stress-frequency 
relationship , the critical shear stress at instability 
is obtained and its variation with radius-thickness ratio 
of the shell is shown graphically. 

Using the technique, similar to that in chapter four, 
the stability of a thin rectangular viscoelastic plate under 
inplane loading is examined in chapter five. The material 
is considered incompressible and the plate is subjected to 
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either axial or biaxial loads. The method of solution is 
similar bo that of Sewell (1973) hut the dynamic effects 
are included in order to study the time-dependent behaviour 
of the material. The critical stress at instability is 
obtained for simply supported edge conditions, When the 
loading is biaxial, the interaction curves are plotted for 
different aspect ratios of the plate. 

The discussion on the results for each investigation 
are placed at the end of the corresponding chapter. However, 
the general observations and all results are summarized in 
chapter six. 



CHAPTER : 1 


CHAPTER : 1 


INTRODUCTION 

1 .1 PRELIMINARY REMARKS 

No mathematical theory can completely describe the 
complex world around us« Every theory is aimed at a certain 
class of phenomena, formulates their essential features, and 
disregards what is of minor importance. The theory meets 
its limits of applicability where a disregarded influence 
becomes important. Thus, rigid-body dynamics describes in 
many cases the motion of actual bodies, but it fails to 
produce more than a few general statements in the case of 
impact, because elastic or inelastic deformation, no matter 
how local or how small, attains a dominating influence. 

For a long time mechanics of deformable bodies has 
been based upon Hooke’s law, i.e, upon the assumption of 
linear elasticity. It is well known that most engineering 
materials like metals, concrete, wood, soil, are not linearly 
elastic or, are so within limits too narrow to cover the 
range of practical interest. Nevertheless, almost all routine 
stress analysis is still based on Hooke's law because of its 
simplicity. 

In the course of time engineers have become increasingly 
conscious of the importance of the inelastic behaviour of 
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many materials, and mathematical formulations have been 
attempted and applied to practical problems. Outstanding 
among them are the theories of Ideally plastic and of 
viscoelastic materials. While plastic behaviour is 
essentially nonlinear (piecewise linear at best), 
viscoelasticity, like elasticity, permits a linear theojry. 

The investigation of stability in linear viscoelastic 
material is the subject of the present thesis. Let us 
first briefly review this vast subject of stability, 

1 .2 BRIEF SUBJECT REVIEW 

In invesbigating the stability of a structure, one 
may ask the following questions: 

i) If the system, which is initially in an equilibrium 
state, is given an arbitrary small disturbance, does 
it remain near the equilibrium state? If so, the 
system is stable, 

ii) Does it remain near the equilibrium state, and in 
addition tend to return to the equilibrium state? 

If it does, the system is asymptotically stable, 

iii ) What bounds must be placed on the magnitude of 

the initial disturbances so that, given a disturbance 
within these bounds, the system will eventually regain 
its original equilibrium state? In other words what 
is the extent of asymptotic stability? 
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In practical sense, an equilibrium configuration of 
a mechanical sysbem is said to be stable if accidental forces, 
shocks, vibrations, eccentricities, imperfections, 
inhomogenaeties, residual stresses or other probable 
irregularities do not cause the system to depart excessively 
or disastrously from that configuration. In a mathematical 
sense, stability is usually interpreted to mean that 
infinitesimal disturbances will cause only infinitesimal 
departure from the given equilibrium configuration. The 
more fact that an assemblage is stable in this refined sense, 
does not necessarily signify that it is safe from engineering 
view]Doint, 

There are two methods of stability analysis, one 
is the static method and the other is the dynamic one. An 
equilibrium state is considered stable by the static criterion 
if no adjacent equilibrium (configuration) exists, or if the 
overall potential energy (of the system) is a relative minimum. 
The dynamic method is more general in the sense that all the 
approaches, based on the static concept, are the special cases 
of this approach when inertia forces are neglected. Since 
this method takes into account the inertia forces in its 
formulation, the mass distribution of the system becomes as 
important as the elastic stiffness of the system. The 
response of the system, therefore, becomes a function of both 
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the space and the time coordinate. The method consisbs of 
the following stipulations and conceptsi an unperturbed 
state whose stability is being studied is specified, A 
perturbation is then applied to the unperturbed state so 
that it is transformed into a perturbed one. Certain 
characteristics, called the norms, are emphasized which 
define the states at any desired time. The change in the 
norms during the transition from tho unpertiirbed state to 
the perturbed one under the influence of perturbation is 
determined. Based on this behaviour, a conclusion may be 
reached regarding the stability of the unperturbed state, 
or its instability. 

In examining the stability of structures undergoing 
large deformations, usually two methods have been adopted. 

The first is the method of small perturbations superimposed 
on finite deformation as available in, for example, Biot 
(1965), A system of governing differential equations of 
motion is obtained J these equations together with rate 
boundary conditions and rate material properties constitute 
the rate problem. The second method makes use of a variational 
technique in which a certain functional vanishes in the absence 
of a stable solution ; non-vanishing of the functional implies 
stability, A detailed account of both the approaches# 
especially for quasi-static problems in continuum mechanics, 
is contained in Hill's work (1959# 1961)* 
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The incremental method of analysis developed by 
Biot in the context of the theory of elasticity could be 
extended to analyse the stability problems of viscous and 
viscoelastic media. In fact, this realization opened an 
entirely new phase in the problems of deformation of the 
earth's crust and tectonic folding of geological structures 
(Biot, 1965), The mathematical foundations of the 
constitutive equations of linear and nonlinear viscoelastic 
materials have been derived by many authors. It would be 
sufficienb to mention only a certain group of papers on 
this suboect, namely those by Rivlin and Erickson (1955), 
Rivlin ( 1955 ), Spencer and Rivlin (1959) , Noll (1955, 1958), 
Coleman and Noll (I96O), The general theory of nonlinear 
viscoelastic materials with memory has been formulated by 
Green and Flivlin (1957, I960) and by Green et al, (1959). 

Some aspects of small deformations superimposed on large 
deformations in materials with fading memory have been 
discussed by Pipkin and Rivlin (1961), Zahorski (1965a) 
dealt with the problem of small additional motion superposed 
on fundamental slow deformation of nonlinear viscoelastic 
media and applied the theory to the problem of buckling 
of an incompressible column of arbitrary dimensions (1965b). 
Dost and Glockner (1982) also investigated the dynamic 
stability of a geometrically perfect three element linear 
viscoelastic column using an approximate numerical technique. 
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In fact, very little work is available on the systematic 
investigation of stability of some basic problems in a finitely 
deformed viscoelastic medium. The reason, perhaps, has been 
the difficulty in handling the complex mathematical formulation 
available. It has not been realized so far that the Hill's 
theory (1 961 ) of instability and bifurcation in elastic and 
elastic-plastic continuum, which has been found so useful 
for application in variety of problems (see e,g, Sewell, 1972) 
can be suitably modified to include even viscoelastic 
materials and that the problem can be made identical to one 
of investigating nonconservative problems in the theory of 
elastic stability by the conventional dynamic method. This 
thesis takes up this challenging work and shows the application 
of the theory in four basic problems (Chapters two to five), 
as summarised in the Synopsis. 

1 ,3 RATE PROBLEM 

In the rate problem the stress distribution, the 
deformabion state and the constitutive material properties 
are supposed to be known at a generic instant and the task 
is to determine the ensuing motion of the body under an 
incremental change of the external loading, 

A common type of boundary conditions for a mixed 
problem is that the rate of change of nominal load vector, 
which is independent of any structural geometry changes, is 
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prescribed on some part of the surface of the body and 
particle velocities on the remaining surface , As the 
current traction is then related to unit reference area, it 
is convenient to introduce the nominal (or, the Lagrangian 
or, the first Piola-Kirchhoff ) stress tensor. 


Referring to a fixed coordinate system, v/hich for 
simplicity is taken as Cartesian here (Figure 1.1), the 
coordinates of a particle orginally at in some reference 
(initial) configuration is denoted by in the current 
configuration B at time t. The vector areas of the same 
surface element are then denoted by n^dS° and in the 

two configurations, respectively. The current load vector 
acting on a surface element is 

dP° = T° dS° ... (1.1 ) 

v;here T? is the nominal traction. 

The nominal stress tensor is then defined by 


ipO = a 


( 1 . 2 ) 


This means that s^^ denotes the current j-'th 
component of the load per unit reference area, acting on a 
surface element which has the normal in the i*“th direction 
in the reference configuration. 
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Now the true traction acting on the surface 
element is given by 

dP^ = as ... (1.3) 

where by Cauchy's stress hypothesis 


n. 


a 


... ( 1 .^) 


0 ■ i iO 

*^ij true (or the Cauchy) stress tensor. Then 

n. 0 dS = n° s, dS° ... (1 .5) 

1 k ka 


The vector areas of the surface element in the two 
configurations are related through the basic formula (cf. 
e.g. StorAkers, 1973) 

Q O Q Y 

n dS = - — , — ^ . n? dS° ,,, (1.6) 

1 P ^ 

where and P are the mass densities in the reference 
and the current configurations respectively. 


Introduction of (1,6) into (1.5) then yields the 
basic relation 


Again 


p° ^ 

^kl " p ax 

d 




... (1.7) 


* « * 


( 1 . 8 ) 
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where is the Kirchhoff stress (or second Piola- 

Kirchhoff stress) tensor. It can be sho’/m that 


a. . 
10 



... (1.9) 


Similarly, one obtains 





... ( 1 . 10 ) 


If the coordinate is chosen to coincide with 
current instant 


X 


i 


such that 


O'. 

10 



P 








at the 


... ( 1 . 11 ) 


Let us now consider the different stress-rates which 

are commonly used. First is the material derivati've, i.e, 

D( ) 

the time derivative following the element, denoted by 
or simply by ( * ) . Again 


D( ) a( ) a( ) 

„ — — — . ,,, ( 1 . 12 ) 

Dt at ^ dXi 

where v^ are the velocity components. 

Other two most commonly used stress -rates are the 
Oldroyd or the convected derivative (Oldroyd, 1950), which 
is associated with the axes fixed in the body and also 
deforming with it, and the Jaumann derivative which is 
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associated with the body and rotating with the body, but 

<5 ( ) 

not deforming with it. The former is denoted by ■ ■ 

3y( ) 

and the later by ^ They are expressed as 

D 


6t ^ ^ ^id ” ^ ^ik ^G,k “ ^ ^jk ^i,k *•' 


d { ) 


- ( )j.j^ ... (1.14) 


where ^ as the antisymmetric part of 

the velocity gradient tensor and a comma signifies partial 
differentiation with respect to Xj_. 

Equations of Motion and Boundary Conditions 

Neglecting body forces the condition for balance of 
momentum is 


/ dS° = 0 


... (1.15) 


in the reference (initial) configuration, or introducing s^^ 

/ s,, n° dS° = 0 ... (1.16) 

gO 1 

Application of Gauss' divergence theorem then yields the field 
equations for equilibrium of the current state 


« 0 


... (1 .17) 
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and for continued balance 


dXjL 


= 0 


... (1.1 a) 


where a dot denotes material rate of change. 


When the inertia forces are considered, the 
corresponding equation of motion in the Incremental form, 
will be 



••• (1.19) 


where Vj are the particle velocities. In fact, v. should 
J ^ 

be interpreted as the incremental displacement so that, at 
the onset, the objection in units of the r.h.s. of (1,19) 
does not arise. 


It follows from the law of balance of moment of 
momentum applied to the current configuration that the true 
(Cauchy) stress tensor is symmetric, which implies for s^^^ 
through (1.7) 


3x. 




no ax. 


n 


... ( 1 . 20 ) 


Material derivation then yields 


Qx. 


ev. 


=mi 


ax 


+ s 


m 


mi 


dx. 


k 


Si 

■ 9X 


m 


. ax. 3v. ex, 

®nj ®nj ^ 

(1 .21 ) 


• • • 
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where are the particle velocities. Material derivation 
may be performed with respecb to any monotonically increasing 
parameters which does not have to be real time except when 
dynamic effects are considered. 


When formulating the rate problem it is convenient 
for the present purpose to leb the reference (initial) 
configuration and the current configuration coincide at the 
generic ir^’tant. The equations of motion (1,19) then 
simplify bo 




« P V 


... ( 1 . 22 ) 


and 


+ s. 


’ij 


V. 1 
i,k 


= s . . + 


'li* ^j,l 


... (1.23) 


where a comma denotes partial differentiation. 


With this choice of reference frame the true and 
nominal stresses coincide at the considered instant although 
from (1 ,7) their material rabes of change differ through 




ij '^k,k 


0 

jk 



... (1.24) 


When bhe nominal surface ti’acbion rate is proscribed 
on Sy and the particle velocities are spenified on , the 
dynamical boundary condibions are 



• • « 


(1,25) 
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In the present context only the instantaneous 
motion is of interest. As a measure of continued deformation, 
the rate of deformation tensor 


e 


i3 



+ V 



.. (1.26) 


is adopted. The vanishing of this tensor is a necessary 
and Sufficient condition for instantaneous local rigidity. 


Constitutive Equations 

Viscoelastic materials show an increasing deformation 
under sustained loads, the rate of strain being dependant on 
the stress. These materials sometimes show a pronounced 
influence on the rate of loading, the strain being larger 
if the stress has grown more slowly to its final value. The 
constitutive equations of these materials may be linear or 
nonlinear. Here we mention only two fundamental models (i,e, 
the Maxwell and the Kelvin models) of linear viscoelastic 
materials. 

The simplest Maxwell model (Figure 1 ,2) consists of 
a linear spring and a dashpot connected in series. The 
extensional strain S of the spring follows the relation 

a = E e ... (1.27) 

cL 

while the extensional strain 6-^^ of the dashpot obeys the law 

ee^ 

o- = — - = pe ... ( 1 . 28 ) 

at b 





where o j_s the load per unit area, E and F are material 

constants , G and 6. are the strains in the spring and 

D 

* 

dashpot respectively, 6^ is the time derivative of the 
strain. Since, both elements are in series, the following 
relation is obtained 

' ' ' ' “ + = 8 + 6 , = 6 

E F a b 

or, alternatively, 

a + qg ^ 6 ... (1 .29) 

where , p 

^2 ” = E . 


The behaviour of this model is similar to that 
of a VISCOUS fluid because the final value of stress tends 
to be zero with time. 

As a second example, the Kelvin model for solids, in 
which the linear spring and dashpot are connected parallel 
(Figure 1.3), is considered. It should be understood that the 
stress Of is not to be distributed on the spring and dashpot 
following the law of levers, but at any time the strain of 
the two elements is the same and the total stress o v^in be 
split into cfg (spring) and (dashpot) in whichever way it 
is necessary to make 6 the same. When applied to this model, 

a = F e 
b 


and 





15 


Further, since a = » 

0 = E e + F e 

or, in the standard form, 

d 

O , (^o + ™ ) e ... (1 .,0) 

where = E and = F , 

The hehaviour of this model is similar to that of 
viscous solid because stress tends to be finite, obher than 
zero, with time. As sqch, normally bhis model is preferred 
to investigate the stability in viscoelastic materials. 

The response of an element of viscoelastic solid to 
an infinitesimal change in the stress is assumed to be of the 
form 

Change in stress = f (Change in strain) 


or, in terms of sbress-rate and strain-rate components, 


Stress-rate component = f (Strain- rate components) 

... ( 1 . 31 ) 

It is desirable that the defining stress-rate be independent 

of any rigid body rotation imposed on the element. The Jauraann 

derivative of the Kirchhoff stress " -w— , which is associated 

t 

with coordinate axes having the same spin as the material 


element and instantaneously coincident with the fixed coordinate 
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axes, is a suitable stress-rate. Consequently, following 
Hill ( 1959 ), ( 1 . 31 ) can be written as 


JA 


D t 


= K. , T e , 
ijkl kl 


... (1.32) 


where the K. , is the tensor of the viscoelastic moduli 
1 jkl 

whose most general form is (Biot, 1965): 

5 ^ 

%Clkl ' ’’ijkl ~ * f ^loki '^1' 

o 


... (1.33) 


subject to the requirement that 


^ijkl ” ^Oikl ” ^ijlk ' ^lij 


(1.34) 


If the incremenbal isotropy is preserved, and if the 
viscoelastic material is assumed to be incompressible and 
the material is supposed to possess the property of fading 
memory = 0 in 1.33), the tensor shown 

to have the following form (Biot, 1965): 


K. . 


ijkl 


(Q. (6 6^^) .X 6„6 


9 1 


ij kl 


(1.35) 


in which Q, ??, A are the material parameters which are 
functions of the current stress; and 6^^ is the Kronecker 
delta. 
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as 


In view of (1.35) , equation (1,32) can be rewritten 



D t 


d 

(2Q + 2r] — — ) e. . + 

at 



6 . . 


( 1 . 36 ) 


In the above, considering incompressibility, i.e,, 


e 


ii 



0 


... (1.37) 


the quantiby \ 6^^ is indeterminate and equal to p, an 
arbitrary' additive hydrostatic pressure. Hence, finally 


D TT . . d 

si = (2Q + 27) — ~ + p . 

Dt at 10 ^ 

in which rj is the viscosity coefficient. 


ho 


. .. ( 1 . 38 ) 


It is worth mentioning here the connection between 
the nominal stress-rate §. . and the Jaumann derivative of 

n 'C . ^ 

the Kirchhoff stress — -LO which is as follows 

St 




+ a 

ik * 





... (1.39) 


where, 6 = — (v. , + v^ < ) is the strain-rate and is the 

10 2 i»0 0»i 

symmetric part of the velocity gradient v. 

J-p J 


The relation (1.39) is extremely important in view 
of the fact that the equations of motion (1,22) are expressed 
in terms of nominal stress-rate For the sake of convenient 



reference in the following chapters, the constitutive equation 
( 1 . 3 s) can now be expressed in terms of s. . as follows: 



(2Q + 2t] 




+ 


a 

ik 



^ f3 , 
;3k ik 


JD'^ . 

It is to be noted that while 

D t 

is not symmetric. 


... (1.40) 

is a symmetric tensor, 


1 .4 STABILITY CRITERION 


Let the velocity v^ be prescribed on some part of 
the surface S of the body and the dead loads be specified on 
the remaining surface S^ , The body which is in equilibrium, 
is imagined to be perturbed further by some kind of external 
agency. The current configuration of the body is said to be 
in a stable state if the magnitude of the ensuing displacement 
decays with time and is vanishingly small, whenever the 
perturbation itself is smal] . If, on the other hand, the 
amplitude remains finite for one type of disturbance, however 
small this might be or tends to increase with time, the state 
is said to be unstable. Therefore, in general, a sufficient 
condition for stability should be that the internal energy 
(the kinetic energy and the internal potential energy) stored 
should exceed the work done by the external loads during the 
resulting motion. 
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In the time interval 6t , due to a virtual displacement 
V. , <5t, the nominal stress changes from s . . to (s.. + s. ,6t), 

•L U u ij 

still referred to the current state. Therefore, the increase 
in the energy will be , 


6t 


/ 

V 


(s . . 

^ 10 


+ 



6t) V . dV 
0>i 


... ( 1 . 41 ) 


The kinetic energy stored is 

/ ^ p v^ dV = 1 / p (6^)2 

The virtual work done by the external forces during 
the time interval 6t, under the load on surface and rigid 
constraints on surface , i,e, v^ = 0 on , is 


6t I T. V. dS = 6t / n,. s. . v. dS = 6t 


S 


10 




^i3 


The total energy stored in the system minus the 
external work is 

(6t)^ [ / 1 s V dV + / 1 P dV ] 

y2 10<J>-^ y2 0 

which should be positive for the system to be stable J hence, 
for stability, 


''j.i ^ J P ] av > 


0 


( 1 . 43 ) 


for all continuous differentiable velocity fields which vanish 
on (but are not identically zero) and are compatible through 
the constitutive equations with the stress rate 
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Fig.1’2 Spring and dashpol in senes 

Maxwell material 



Fig,1'3 Spring and dashpot parallel 

Kelvin material 

I 
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CHAPTER : 2 


INSTABILITIES IN RECTANGULAR VISCOELASTIC 
SOLIDS UNDER PLANE STRAIN 

2.1 PROBLEM STATEMENT AND ITS FORMULATION 

An incompressible viscoelastic rectangular solid 
undergoing continued deformation in plane strain is 
considered. The dimensions in the initial configuration 
are 2a^x 2b^ x 20^, It is brought to the current configuration 
by means of finite deformation (bension or compression) along 
one direction such that the dimension 20 ^ of the specimen 
remains constant, i.e, the specimen undergoes deformation 
under plane strain. In the current configuration, the 
dimensions of the specimen are 2a x 2b x 2c^. Referred to 
this state, a fixed Cartesian coordinate system x^^, coinciding 
with the axes of symmetry of the specimen, is chosen as a 
reference frame. Whenever convenient, the coordinate x^ 
will be replaced by x,y,z and, the velocity components by 
u,v,w, respectively. With respect to this frame, the 
deformation is assumed to be parallel to x-y plane, and the 
distribution of stress is assumed to be 

c. - = = s.. = 

iO iG 




a 

0 

0 


0 

0 

0 


0 

0 

a 

2 


( 2 . 1 ) 
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Under this state of stress the current configuration 
is in equilibrium. The aim is to find out under what 
conditions the current state is stable/unstable. 

During perturbation from the current state, it 
is assumed that the longitudinal ends x = ^ a are frictionless; 
hence the shear fraction rates on these ends will be Eero, 

Also, the lateral faces y = + b are supposed to be free of 
nominal traction rates. These conditions are expressed as, 

= ®i 2 = 0 , X « jt. a 

... ( 2 , 2 ) 

= ®ij = 0 > y = ± b 



where n^ is the unit outward normal to the boundary surface. 

For the present case, the rate consbitutive equation 
(1 ,40 ) yields following relations: 

d 


'11 

(2C1 + 2i7 — 

ut 

-O) + P 



22 "" 

(2Q + 2?7 — ) 
dt 

■^2,2 P 



’12 " 

a 0 

" r 

) ^2 -v (Q + 

0 

7, + 

■Ot 

— .-) V, 
2 

21 

a 

(Q + 77 — ■ 

at 

7 ) Vi ^2 + 

a 

+ 77 -Sf"" 

' Qb 

a 

2~ 




2,1 


> 


V, 


2,1 


> » « 


(2.3) J 
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The equation of motion (1.22) now reduces to 


11 ,1 


= P 

« ^ 

^1 




• f 

CM 

®22,2 

- P 

^2 


... (2.4) 


Finally, the incompressibility condition (1 .37) 
for a plane strain problem may be written as 


^1,1 ^ ^2,2 = ° 


... (2.5) 


2.2 INSTABILITY COTOTIONS 

Substitution of relations (2.3) in the equation of 
motion (2.4) yields the following two equations in the 
physical components of the velocity u,v and the parameter p 

0 0. o 


+ (2Q +2 77 


0. a 

- a) u + (Q + T) ) u 

^ XX ^ ' 5 1 2 yy 


0 a ^ a 

p + (Q + 7] — - -r ) u -h (Q + 77 -r- + “ ) V 

y ot 2 3t 2 ^ 


... (2.6a) 


+ (2Q + 2 77 — — • )v = P V 

ot yy 


... (2.6b) 


where, now a subscript denotes the partial differentiation. 
The boundary conditions (2.2) at y = + b yield 5^2 = C) 


and =0 i.e. 


(2Q + 2r) ) V + p = 0 ’'I 

et y 1, 


u + V 

y x 


y = + b 


... (2.7) 
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Let the velocity field, and the parameter p be 
assumed in the following form 

u = E f(y) sin kx exp (wt) 

V = E <P(y) cos kx exp (wt) ... (2.8) 

p = £ ■'g(y) cos kx exp (wt) 

where k = rox/a and w = i w i = V --I , u being the 

frequency, f(y), ^P(y) and g(y) are functions of y to be 
determined. The assumed velocity distribution satisfies 
boundary conditions (2.2) at x = + a. Substitution of 
equations (2.8) in the incompressibility condition (2.5) 
and the equation of motion (2.6) yields three ordinary 
differential equations in f(y), d(y) and g(y) : 

= 0 ... (2.9a) 

kg - (2Q h 2p w - a + ~ ) + (Q + nw - ^ ) f " 

-(Q + 77 u - I ) k(p’ = 0 ... (2.9b) 

, /n - P 9 T, 

“g -(Q + 77W+« + — )k^<P-h(2QH2rjm)<P 

2 k2 

+ (Q + ^ 5 - I ) kf’ = 0 ... (2.9c) 

where, now, the prime denotes differentiation with respect 
to y. 
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Eliminating f and g from the set of equations (2.9), 

the following fourth order equation in <P{y) is obtained: 

-2 

- a lilt _ p U) 

(Q + r)(0---)<P “ (2Q + 2rf (ji + ) (p' ' 

2 1^2 

-2 

a P ^ / 

+ (Q+ t] (li + — + ~ \ <|) = 0 

2 k2 

which, alternatively, can be written as 


(1 +\co “ ©)'^ - (2 + 2X td + 


) k^<p' 


-2 

L 

+ (1 +\a) + 0 + ) k^ tp = 0 

k2 


... ( 2 . 10 ) 


in which 0= — JL ^ a 

2^ Q Q 

With the help of equations (2,8), the boundary 
conditions (2.7) may' be written as 


(2Q + 2?7a))V ~g = 0 


f ’ - k<P 


= 0 


y 


... (2.11) 


which after eliminating the functions f and g, take the 
following forms 

r , o'" 

(3 + 3 - 0 + — “ )<P - (1 +X w - 0) = 0 ... (2,11a) 


» « ,2 
(p + k (P 


... {2,11b) 


such that (1 +X (3 - 0) ^ 0 , 
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Referring to the equation (2,10), let the solution 
of <P be taken in the following form 

(p (y) =s A exp C ^ ky ) 

where k = nn/ a as before. 


... (2.12) 


The characteristic equation is 


r^ - rp = 0 

in which 


... (2.13) 


r^ = (1 + X m - 0 ) 

-2 

am 

Tp = (2 + 2\ m + — — ) 

k2 

a 

p = (i+Xm + 0 + ) 

^ k2 

The roots , '^2 » and of the characteristic eq.(2,13) 


are 


2 V2 

_ 3^2 i ^^2 " ^ ^ 1 ^ 3 ^' -*^/2 

i9i = i [ ^ ^ 


] 


2r. 


1 




and are ordered such that ^2 “ " ^4 

On Simplification, 


... (2.14) 


-2 
a oi 


^ ~ B 1 , and 

W 2 - - ^4= [ T-*- ^ ] 

1 +XW - 0 


... (2.15) 
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Two distinct cases arise, viz,, 

i) Case of equal roots i.e, when ^ = ^2 
ii) Case of unequal roots i.e, when p ^ 2 * 

2 , 2,1 Case of Equal Roots 

Equal roots are possible only if ^ = - ^ 3 =“ “I* 

i.e, when ^3 expression of in equation (2.14), 

This leads to 

a 

0 + ^ Q 

2 k^ 

The solution of V (y) for this case is given by 

<P(y) a cosh ky + sinh ky + sinh ky 

+ A^y cosh ky ( 2 . 16 ) 

where A^ ( i = 1 , 2 , 3 , 4 ) are arbitary constants. 

Substitution of the boundary conditions ( 2 , 11 a and 
2 , 11 b) at y = + b leads to four homogeneous, linear, algebraic 
equations for the constants A^ ( i = 1,2, 3, 4). For a non- 
trivial solution to exist, the determinant of the co-efficients 
A^ should vanish J thi t leads to 

a ^1 a ^2 ®13 ^14 

®11 ^12 “^13 “^14 

^21 '^22 ®23 ®24 

^21 “^22 ^23 ^24 




0 


(2.17) 
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where 


^ =2 cosh kb 
2 

a ^2 " cosh kb + 2 b sinh kb 

a,|^ = 2 sinh kb 
2 

a. i, = — sinh kb + 2 b cosh kb 
^21 " ^ sinh kb 

^22 “ "■3C2) sinh kb + (c^ ^ ^ 2 ^ cosh kb 

®^23 “ "■ ^ 2 ^ ^ cosh kb 

^24 “ ^*^1 " ^°2^ cosh kb + (c^ - C2) kb sinh kb 


in which 

C2 


3 + 3 \ (o-0 + 
1 + \ io - 0 



k2 


The simplification of the condition ' 2 . 17 ) results in 
the following relation 


(2kb) 

■ ■ I 11 ■ Ml ■ ^iii 

sinh (2kb) 


... ( 2 . 18 ) 


It can be easily shown that the condition ( 2 , 18 ) 

b 

is satisfied if, and only if, “'^0 • Hence, instability 
is not a possibility in case of equal roots. 
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2,2,2 Case of Unequal Roots 

In case of unequal roots (i.e, '3^^ / the roots 

of the characteristic equation (2,13) "being obtained from 
the expressions (2,14 or 2.15), the general solution of ^(y) 
may be written as, 


‘^^(y) = cosh 3^ ky + cosh ^2 ky + 
+ sinh ^2 ^^7 

where, A^ ( i- 1,2,3»4) are arbitary constants. 
Substituting, 3-j = 1 , and > 

“2 
a 0 ) 

1 + A 5 + e + — ^ 1/2 

3 = “^2 " [ — ::: — ^ — ~ ^ 

1 + X tS - 0 

equation (2.19) is rewritten as 


sinh ky 
... (2.19) 


... ( 2 . 20 ) 


V(y) B A^ cosh ky + ^2 cosh 3 ky + sinh ky 

+ A^ sinh 3 ky ... (2,21) 

Substitution of the boundary conditions (2,11a and 
2,11b) at y = + b leads to four homogeneous, linear, algebraic 
equations for the constants A^ ( i = 1,2,394). For a non- 
trivial solution to exist, the determinant of the coefficients 
should vanish as before. The condition for instability 
may thus be put in the following form'. 
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a^1 


^2 

^13 

®14 

"^11 


“^12 

®13 

^14 

^21 


^22 

^23 

^24 

^21 


^22 

"^23 

'^24 

noWj, 

^11 " 


Sinh 

kb ; 

^21 " 

^12 " 

‘^2 

sinh 

/3 kb; 

^22 " 

a^3 == 

^1 

cosh 

kb ; 

®23 

^14 

^2 

cosh 

/3 kb; 

^24 


= 0 


c, cosh kh 
D 


= cosh ^ kb 


0 ^ sinh kb 
c^ sinh /3 kb 


( 2 . 22 ) 


in which 


a 


= (3 + 3kW“© + - — “) - (1 +\ w - e) 

k2 

-2 

, ^ ^ aw - 

02 = (3 + 3 ^ “ 0 + )j3 - (1 +X w - 9) 


“ 2 , = ( 1 + ^ ) 4 


The 4x4 determinant (2,22) can be decomposed into 
the product of two 2x2 determinants, viz., 

A-j .4^ « 0 ... (2.23) 
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where 


and 


11 

a^2 

^21 

^22 

'13 

^14 

'23 

^24 


(2.24a) 


(2.24b) 


Hence, for a non-trivial solution of A^, equation (2,22) 
is satisfied if , 


either , 

Ao 


or 


0 

0 


• * • (2 ,25a) 

... (2,25b) 


The above decomposition has a simple physical 
inberpretation J the conditions (2,25a) and (2,25b) correspond 
to instability phenomenon in the antisymmetric and symmetric 
modes, respectively. 

For instability in the symmetric mode, ^2 = 0 (cr.2.25b) 


leads to 


tanh (j8 kb) 


tanh kb 


4i0 

(1^^^)^ 


(2.26a) 


For instability in the antisymmetric mode, 
4| = 0 (cf. 2,25a) reduces to 


tanh (|3 kb) ''' •-2\2 


tanh kb 

where ^ is given by (2,20), 


(2.26b) 



33 


2.3 BEHAVIOUR OF THE INSTABILITY CONDITIONS 

In addition to 0, (the non-dirnensional stress 
o 

parameter — ) ^ the following two non-dimensional parameters 
are introduced? 

„ r , 1 ,1/2 

0 = L “ ] “ ; 2 ] r, ... (2.27) 

^ P Q 


These two parameters are referred to as non-dimensional 
frequency and non-dimensional viscosity in the subsequent 
discussions. Therefore, the instability conditions (2.26a 
and 2,26b), in terms of these parameters, will read as 


tanh /3 kb 
tanh kb 

and , 

tanh (3 kb 
tanh kb 


4/3 



(1 + ^ 2)2 

( 1 + 

4/3 


... (2.28a) 


... (2.28b) 


where, now 2 

Q‘^ 

1 +i^.-Q + 0 — 1/2 

3 =[ -_] 

1 + 0 - 0 

Numerical solution of the equations (2,28a and 2.28b) 
for a particular value of kb (i.e, for a given aspect ratio 
b/a) , will yield various values of /3 (cf. Figs, 2,1, 2.2, 

2.3, 2.4, 2,5), Let the numerical value of 3 obtained by 
solving either of the equations (2,28a) or (2,28b), for a 
particular value of kb, be designated by 3, The value of ^ 



54 


will either be real or imaginary. Actual numerical 
computations did not yield any complex value of 3 ♦ Again, 
it is evident from the figures (Figs. 2.1, 2.2, 2.3, 2.4, 2,5) 


that 






(i) 

if 

3 

is real 

0 < t3| ^ 1 

and 

(ii) 

if 

3 

is imaginary 

i < I pI i, '>'1 



where 7 is real number and is > 1 , 


Once, a particular value of 0 has been determined 
(for a particular aspect ratio), load (stress) vs, frequency 
relationship may be studied from the relation 


3 



1 +i\^ Q + 0 

1 Q 



... (2,29) 


In general, the non'* dimensional frequency n is 
complex J 3 which satisfies either of the relations (2,28a 
and 2.28b) may be real or imaginary J non-dimensional load 0 
is real, either positive (tensile) oc negative (compressive). 
Equation (2,29) on simplification yields the following- 
equation in non-dimensional frequency 0 • 


( - [i(i - ^ ^ [(1-e) (1~3^)+ 20] 0 

n% ‘nit 


from which, 
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Q 


r . -2v n r -P \^n7t ,2 

1 [ (1-13“^) -J— ] i [ flCl-IB^) (-1— ) ) 

n TC '=■ 2 

p 

+ {(1-0) + 20}] ... (2.30) 


The following observations can now be made: 

i) As /3 is real or imaginary (but never complex), 

the second term following the + sign within the 
square bracket in equation (2.30) is a real number, 
positive or negative. 

Further, if the second term within the square 
bracket is positive it represents the real part 
of the non- dimensional frequency. 

ii") The first term within the square bracket in 
equation (2.30) is always positive, since. 

Oil I 1 > when 0 is real and i < j ^ '^i? 

when 0 is imaginary, “y is a real number > 1 . 

iii) When the second term within the square bracket 

is real as well as positive, the first term within 
the square bracket in equation (2.30) will represent 
the imaginary part of the non-dimensional frequency. 
Hence, if the second expression within the square 
bracket is positive, this will lead to oscillatory 
convergence of the system, the non-dimensional 
frequency being complex. 
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iv) 


v) 


form 


If the second term within the square "bracket is 
negative (i.e, the second term itself is imaginary), 
but the second expression itself is less than the 
first term (which is positive imaginary) , the non- 
dimensional frequency in equation (2,30) will he 
imaginary but positive. This will lead to exponential 
convergence of the system. 

Hence, for instability (i,e, divergence), the second 
expression within the square bracket in equation (2,30) 
should be negative, as well as, 

magnitude of the second expression 

2 magnitude of the first expression 
' in equation (2,30). 

Equation (2,30) is rearranged in the following 


= 1 [ (1-ih ( , 


7:2 s 


) ] 


+ i [ {(1-3^) ( tfil) - ( (i-s) (1-3^) + 

2 

... (2.31) 


Hence, the fco-lowing should be satisfied for 
divergence, 

Q iz. 


^ A^nji 2 o , ^1/2 

[ {(1-r) ( ™ ) > -C(i-e) d-r) H 20} ] 


2 




(2,32) 



37 


which on simplification, reduces to 

^ - i ) ... (2.33) 

js^-l 

or ® i. ( ) ••• (2.33a) 

3^+1 

From the requirement (2.33), it is evident that 
for instability jn compression (i,e, 6 negative), in either 
symmetric or antisymmetric mode, modulus of 5 (whether j3 is 
real or imaginary) should lie between zero and one. Further, 
the instability in tension is not a possibility, 

2,3.1 Instability in Compression 

25 56 (l) Symmetric Mode ■ Real Roots : 

Referring to the equation (2.28a) and Fig, 2,1, 

3=0 or +1 always for any value of kb. If kb is 
sufficiently large (i,e, if b/a is large enough so that 
kb > 4 ), another real value of 3 is possible, the absolute 
value of which lies between zero and 0,29561 (cf.Fig, 2.2), 

The value 0 = + 1 need not be considered, since, as has 
already been explained , there can be no instability in the 
case of equal roots. The case 3=0 leads to (cf, equation 2.33) 

0^-1. 0 , 0_£kb^ 4.0 

When kb > 4, the root^ absolute value of which lies 
between 0 and 0,29561 (cf. Fig, 2. 2) will yield a lower value 
of 0 than above. 
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When “■*►00 j i.e, kb — « , equation (2,28a) 

3 . 

reduces to 

Therefore, 0 = + 0.2956'1 and hence 
0 < - 0,83927 
when b/a — *>■ 00 , 

b ' , ^ 

Further, as ^0 , i,e, kb-^K- 0, equation (2.28a) 

B 

yields 

4/3 

(1+|3^)^ 

Solving, 3=0 and + 1 , are the only possible real roots, 
which yield 0 = - 1,0. 

Therefore, for instability in the symmetric mode 
under compression 

- 1 ,0 < 0 < -- 0.83927, when 4 < kb <“ , and 

0 ="1.0, when 0^ kb^ 4.0 ... (2.34) 


Case (li) Symmetric Mode: Imaginary Roots: 


VJhen P is imaginary (let 0 
equation (2,28a) reduces to 
tan(r kb) 4')' 

tanh kb 


= T'i, where 7 is real), 


... (2,35) 
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Solution of the above equation yields infinite 
number of imaginary roots (ref. Fig* 2,4) and equation 
(2*33), taking p - ri takes the following form 

-2 

1 +7^^ 

0 < ^ ( __ ) ( 2 . 36 ) 

1-7"^ 

It is evident that the least value of 0 is -1,0 when 
7=0, which is always a solution of equation (2,35) regardless 
of the value of kb or b/a. Also, for 0 to be negative, 

0 < | 7 | < 1 . 0 . 

Considering both the real and imaginary roots of 
j3 , modulus of which lies between zero and one, the value of 0 
may be stated to lie as follows 

p=0 er=-1.0 1 0<kb< 4.0 

0 < l^j ^ 0.29561, -1 .0 ^.0 <-0.83927; Ai^b < ... (2,37) 

Case (iii) Antisymmetric Modes Real Root * 

Referring to the equation (2,28b) and Fig, 2,3, 

3 = + 1 and another pair of roots j3 , such that 
0,29561 < jisj < 1 ,0,when 0 ^ kb < «» , As before^3 * +, 1 are 
not considered because instability cannot occur when the 
roots are equal, Thi s, for any real value of 3, 0<. j^j^l.O, 

-0.83927 <®<0; 0<kb<~ ... (2.38) 
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As ~-^0 i.e, kb-*-0 and hence, the equation (2,28b) 
reduces to 

(1+3^)^ 

3 = 

43 

whose solutions are 
3 = + 1 

Again as — -*oo ^ i.e, kb -voo , and the equation (2,28b) 

becomes 

^ ( 1 + 3 ^)^ 

43 

whose solutions are 

3 = +, 0.29561, and hence from equation (2.33) 

0 = - 0,83927 

Case (iv) Antisymmetric Mode; Imaginary Roots ; 

When 3 is Imaginary (let ^ ~ y±f y is real as 
before) , equation (2,28b) reduces to 

tan('y kb) (l-r^)^ 

— = - ... (2.39) 

tan h kb 4 y 


Solution of the above equation yields infinite 
number of imaginary roots (ref. Fig. 2,5), but always 
1^2 "1 i*®* iojf’ practical purpose, no imaginary 


roots are to be considered in order to determine the 



41 


non-dimensional compressive load. 


As -■ 


0 i.e, kb -♦ 0 


the equation (2.39) 


takes the following form 


r = ~ 


4r 


... (2.40) 


which does not yield any real value of Y, 


As i.e, kb , the equation (2,39) reduces to 


tan (y kb) = - 


(1-Y^)^ 




... (2.41) 


Solution of this equation (2.41 ) yields infinite 
number of roots modulus of which are greater than 1 , hence 
need not be considered for the compression case. 


Therefore considering expressions (2.38, 2,40, 2.41) 
instability in compression can occur ^ if 


0.29561 <1^1 < 1.0 ; - 0.83927 < 9 <0 ; 0 <kb < “ 

... (2.42) 

Comparing expressions (2,37 and 2.42), it is observed 
that the instability in the antisymmetric mode is the only 
possibility. 
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2.4 NUMERICAL RESULTS AND DISCUSSIONS 

The aim is to find out the critical value of 3 
(real and/or imaginary) from the instability conditions 
(2.28a and2 .28b) and hence determine from the equation (2.33) 
the non‘'dimensional load © at which the material becomes 
unstable. The computations were performed for all possible 
dimensions of b/a. Three values of viscosity (O. 01 , 0.1, 

0,5) were taken. The observations are as follows; 

i) The critical value of the load © is independent of 
the magnitude of non-dimensional viscosity, and the 
nature of motion is divergent. So long as the load 
does not exceed the limit prescribed by equation (2, 33) 
nature of motion will be either oscillatory or 
exponential. As soon as the load exceeds the limit 
prescribed by equation (2.33), motion diverges 
exponentially. Right at the critical value of the 
load, at least one value of the frequency q becomes 
zero i.e, real as well as imaginary part of the 
frequency becomes zero. Although the cricical load 
at instability is independent of non-dimensional 
viscosity the rate at which the motion decays, 

when the load does not exceed the limit specified 
by equation (2,33) > as well as the rate at which the 
motion diverges, when the load exceeds the limit 
specified by equation (2,33), depend on the magnitude 
of the viscosity 4 ^ , This has been verified numerically. 



43 


ii) For a rectangular viscoelastic solid, no loss of 
stability is possible in tension, 

iii) The primary mode of instability in a rectangular 
viscoelastic solid under compression in plain 
strain, is antisymmetric. This means that the 
condition (2,2Sb) yields the lower value of 
compressive load/stress (at which the motion becomes 
unstable) than that obtained from the condition 
(2,28a), This is also evident from the plot of 
nondimensional load/ stress Q versus specimen 
dimensions (b/a or a/b) in Fig. 2.6. 

Further, when “is large (^ = 0 ) , both symmetric 
and antisymmetric modes of instability are possible under 
compression giving identical values of load/ stress 
(e = ^ 0.83927). 




g2A Distribution of real roots in symmetric mode (0 ^ Kb<4) 
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Ffg.2*3 Distribution of real roots in antisymmetric mode 
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Fig. 2-5 Distribution of imaginary roots in antisymmetric mode 


Symmetric mode 



Fig. 2-6 Variation of critical stress in compression with aspect 
ratio in symmetric and antisymmetric modes 
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AXISYWETRIC INSTABILITIES IN VISCOELASTIC 
CIRCULAR CYLINDER UNDER AXIAL LOAD 

3.1 PROBLEM STATEMENT AND ITS FORMULATION 

In this chapter, the stability of a viscoelastic 
solid circular cylinder under axial load has been examined. 

At the current instant, during the process of continuing 
deformation, the cylinder has radius a and length 2L, To 
investigate the instability phenomenon, the behaviour of 
the cylinder in transition from the current state to a 
neighbouring state, under infinetesimal increments of 
boundary values, has been isolated for study. Referred to 
the current state, a fixed cylindrical coordinate system , 
with one of the axes coinciding with the axis of symmetry 
of the cylinder has been taken as the reference frame. 

Whenever convenient, the co-ordinates x^, Xp,x^ will be 
replaced by r,0,z, respectively, and the velocity components 
v^>'V 2»^3 respectively. 

In the current state, the internal stress distribution, 
assumed homogeneous, is supposed to be given by 

0 0 0 

a. . = 0 0 0 

0 0 ° 
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where cr is the stress in the axial direction. 

During the incremental deformation from the 
current state of equilibrium, the lateral cylindrical 
surfaces are assumed to be free of any tract ion*- rates 
while the cylinder ends are frictionless (i.e, shear 
traction-rates are zero on these two ends), and these two 
ends move with constant velocity relative to each other. 
Therefore, the boundary conditions for the non-homogeneous 
incremental deformation are 





^3 

• 

T . 

D 

= 0 

= ’^i 

• 

s. 

13 

••• 

- 0 , a = 1 ,2 

(3.2) 


• 

T • 


"i3 

= 0 , r = a 

(3.3) 

and, 

• 

T 

^3 

= 

®i3 

are finite at r = 0 ... 

(3.4) 


where Is the unit outward normal to the boundary surface, 

Alternatively, these boundary conditions may be expressed 


as follows 



w = 

0 


2 - + L ► 

... (3.2a) 


# 

s = 
zr 

s^e 

= 0 

* 



■ 

®rr 

Sre 

• 

‘-’rz 

= 0 , r = a 

... (5.3a) 

and , 

• 

®rr » 

Sr0 » 


are finite at r = 0 

... (3.4a) 


.-?^.NTRAL library 


^ 
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For axisymmetric case, u » u(r,E) , v = 0, w = w(r,z), 
and the equation of motion (1 ,22) assumes the follov/ing form 
in the present case 


• I 

®rr,r ^ ®rr ~ ®2r,z ® 


• 1 • 

®rz,r " ®rz 


(3.5) 


= Pw 


2^ j. ^ ZZ ^ Z 

where, as before, p is the current density of the material, 
and a comma indicates partial differentiation* 


Furthermore, the condition of incompressibility 


(I.S?) reduces to 

u ^ + w _ =0 (3.6) 

» ^ r » 


For the axisymmetric case considered here, the 

1 

components of the strain-rate tensor 6=— (y_ +v. .) 

2 1 , j j , 1 

•1 

and those of the rotation tensor m . = r- (v. , ~ v. ) are 

ij 2 ' 1,3 3,1 


given as follows 


e 


id 


h 

ar 


0 

1 

2 


* 

a z 


^ ) 

ar 


0 

u 

r 

0 


2 ez dr ^ 

0 

dw 

a z 


P • i 


(3.7) 
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0 


CO 


la 



9u aw 
d z 9 r 


0 

0 

0 


1 ^ ^ Qw ^ 

2 92 5r 

0 

0 


With the help of the relations (3.7) and (3.8), 
constitutive equation (1 ,40) yields following relations 


= (2Q + 2riD) + p 

Sgg = (2Q + 2r) D) - + p 

r 

®zz + 2?7D-a)w_+p 

s = (Q + 71 D “ “• ) (i) + w ) 

rz ' ' 2 ^ ' ,z ,r^ 

®zr = (Q + 7 ? D + ~ ) + (Q + 77 D - I ) w^^ 

a 

where D stands for the operator , 

at 


\ 


(3.8) 


the 


(3.9) 
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3.2 INSTABILITY CONDITIONS 

Substitution of the relations (3.9) in the equation 
of motion (3.5) results in the following two equations: 




(2Q + 2 77 D) u^^^ + P,r ^ (20 + 277 D) (u ^ ~ ) 


+ (Q+r7D+^)u,^ + (Q+77D-“)w^2 = Pu 




and, 


Q* ^ CJ 

(Q + 77 D “ - )(u + w ^^) + - (Q + 7) D - “ ) (u ^ 


,rr 


+ (20 + 277 D - n) w __ + p 


zz 


J 


= P w 
... (3.10) 

In view of the boundary conditions (3.2a), the solution 
of equations (3.10) may be taken in the following form 


u = f^(r) cos kx exp ( wt ) 
w = f2(r) sin kx exp ( wt ) 


... (3*11) 


p = f,(r) cos kx exp ( cSt ) 

J 

n - 

in which k = •“ — , n being an integer, (o = i«, 

L 

w being the frequency, i = -1 , and f^(r) are unknown 

functions to be determined. 
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Use of (3.11) in (3.10) and in the incompressihility 
condition (3.6) yields the following three equations in 
f ^ , f 2 and f ^ : 

■PI - o P P 1 

3; - (Q+ 770) + -+ — ) (2Q + 277 w ) f. 

+ i (2Q + 2 77 6 ) ) f^' + (2Q + 2 77 u) f^' '+ (□ + 77 5 - ^ )kf 2 

= 0 ... (3.12a) 

- kf^ - (Q + r7 oi ~ I ) kf^’ - ~ (Q + r7 to - | ) kf^ 

f r\ — CJ . t t 1 ^ rr I 

+ (Q+77(0--)f2 +-(Q+77W--)f2 

-2 

9 ^ 2 

- (2Q + 2 77 to - a + — ) k f^ *= 0 ... (3.12h) 

k2 

t ^'\ 

~ + kf 2 =0 ... (3.12c) 

in which, now, the prime denotes differentiation w.r.t, r. 
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Eliminating ±2 ^5 above equations, the resulting 

equation in f-|(r) is arranged as follows: 


T f I 1 


f t f 




1 V 


+ a2 + a, + a,, = 0 


^3 "1 


4 


... (3.13) 


whe re , 


a. 


a. 


L 

r 


-2 

aw 


(1 + xw + —ly ') 

2k 2 [ ii- ] . 3. 

(1 + XW - 0) 

^2 

a w 

2 k^- (1 + xw + — -2 ) 

— - ■ ^ ^ 2 Ic 


(1 + XU - e ) 

-2 

a 0) 

(1 + XU + 0 + — — ) 
[ 


] + 


(1 + XU - 0 ) 


2k‘ 


(1+ Xu + 


a u^ 


j- _ 2k^ 

(1 + XU - 0 ) 


) 


L, 

r3 


3 

r4 


such that (I+XU-0) 0 , 

Again, as in Chapter 2, 


Let, 


_“a32 ° 

(1 + ^ ) 

t) = ~ — — 


(i + \0 - 0 ) 


(1 + +0 4- 


-2 

a 03 


c = 


(1 + \03 - 0 ) 
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Observing that (1+c) =~2b, the preceeding equation 

(3.13) nay be written as 




! 1 t f 2 ? » ? 

r ^ 


t I 


r2 

1 


^1 ^ 


- 


3 


+ 2k2b [ -* — ] + k%f^ =0 ... (3Jk) 


t f 1 

r “ r2 
Introducing the operator 


l2 = 


d 1 d 

C “ * T” (i'...) ] = ( 


dr 


r dr 


the equation (3.14) simplifies to 


d^ 1 


dr^ r dr 


(L^ + 2k^ b c ) f. 


or 


1 


(L^ + (L^ +k^/32) f^ 


- 0 

= 0 


... (3.15) 


in which 3^ and ^2 represented a£ 


3^ , 02 = [ b + Y (h^^c) ] 


1/2 


Noting the relation 2 b-- (1 + c), it is evident that 


0^ = Y-* c , ^2 " ^ ••• (3.16) 

The assumed velocity field (3*11) identically satisfies 
boundary conditions (3.2a). With the substitution of this 
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velocity field and the use of the constitutive relation 
( 3 . 9 )# the boundary conditions (3.3a) become 


( 2Q + 2r) oi ) + f:^ = 0 

( Q + 77 5 - I ) ( - kf^ + f^) = 0 


r - a 


... (3.17) 


which, after making use of the equations (3.12b and 3.12c) 
assume the following form 

do p df/i n 

— (r L'^f.) + ( 2b - 1) k^r — S. _ , 0 

dr ^ dr ^ 

( ) f^ = 0 



... (3.18) 

in which, ^ -2 



( -] + - 0) 

and, b and c have already been defined. 

Two distinct cases may arise in the solution of 
the equation (3.15) J 

(i) P 2 and (ii) ^ • 
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Case 1 i ^2 

The general solution of the equation (5 #15) may be 
written in the following form 

f^ *= ( k /3^r) + ApJ^ (k ^2^) 

+ ( k |32r) ... (3.19) 

where, ( i = 1 >2,3,4) are the constants of integration 
and, J-j and are the Bessel functions of the first and 
second kind , respectively. Since f<^(r) should be finite 
at r = 0 , A^ and A^ must vanish, and equation (3.19) 
reduces to 

f^ = ilj ( k /3^ r) + A 2 J^ ( k ^ 2 ^) ••• (3.20) 

The constants A^ and A 2 are to be determined from 
the boundary conditions (3.18), Thus, by substitution, 
the following two equations are obtained? 

[ - (2b + e) ( k i3^a) - k B^a (l-2b + i3^ ) ( k /3^a)] A^ 

+ [ - ( 2b + e) (k ^ 2 ^) - k 32^ ( 1“2b + ^l) ( k 32^)1 ^2 

= 0 

[ (1 - 3f ) ( k 3^a)] A, + [ (1 - 32 ) J>i ( k32a)] 42 = *^ 


• * ^ 


( 3 , 21 ) 
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For instability, there must exist a non-trivial 
solution of (3.2'1) , i,e, the characteristic determinant 
of the co-efficients and must be equal to zero. 

Since 02 = i and 0^ = ^[-c and, in general should be 
taken as a complex quantity, the characteristic equation is 

(1 + c) ( k 0^a) (2b+e) (k 0^a) + 2k 0^a J^(k 0^a) 

^ ^ ^2^^ (2b+e) (k 02 ^^ (1 +c) k Pga J^(k 02 a) 

... ( 3 . 22 ) 

in which, b, c and e have already been defined earlier. 

Case II 5 0^ = 02 = i* 

The general solution of the equation (3.15), non- 
singular at r = 0, may be written as, 

f^ = Ajl^(kr) + A2r ^^(kr) ... (3.23) 

where, I^ and I are modified Bessel Functions of first kind, 

J o 

A similar procedure as in the Case I, leads to the 
following condition for instability: 

[ 1 + (kr)^] I^(kr) = (kr)^ iQ(kr) ... (3.24) 

The above condition is satisfied for kr = 0 only 
(i.e, r/L = 0), That it is so can also be seen graphically 
as in Fig. 3.1. Hence, the loss of stability is precluded 
whenever 0^ = 02 * 
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3.3 NUMERICAL COMPUTATIONS 

As in Chapter 2, in addition to S (the non-dimensional 
a 

stress parameter ), two additional non-dimensional 

parameters are introduced to facilitate numerical computations. 
These are 

r P 1/2 1 1/2 

Q - L ] 05 ; \ = [ J rt ... (3.25) 

Q ^ p Q l2 

These two parameters are referred to as non-dimensional 
frequency and non-dimensional viscosity in the subsequent 
discussions. Hence, in equation (3,16) will read as 

2 

1 + i X. Q + e - 

[ ^ ] ... (3.26) 

1 + i X^ fi - 0 

The equation (3.22) expresses the relation between., 
the non-dimensional stress © and the frequency parameter q 
which , in general, is a complex quantity. It should be 
noted that , for a particular value of the stress, if the 
imaginary part of the frequency Q is positive, the motion 
is stable and, if the imaginary part is negative, the motion 
of the system becomes unstable at that stress J this is 
regardless of the sign and the magnitude of the real part 
of the parameter Q « Hence, it may be stated that the value 
of stress will be termed critical below and above which the 
imaginary part of the frequency parameter is positive and 
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negative, respectively. It has been observed that, during 
this transition, the real part of the frequency is very 
small. Such a calculation is sho-wn in Table 3.1 for 
a/L = 1 ,0 J = 0,5, For such a case, the critical value 
of 6 comes out to be - 2,919896, 

TABLE 3,1 : VARIATION OF FREQUENCY WITH STRESS NEAR THE 


CRITICAL STRESS FOR X.=0.5; 0^ = -2.919896 

j c r 


Stress Parameter 

0 

Frequency Parameter 

Q 

Remarks 

- 2.819696 

1,21x10“^, 6.92x10''^ 

Below critical 

- 2.919696 

3.35x10"''°, 6.59x10"^ 

Belov/ critical 

- 2.919896 

15 n 

1.68x10“ , 2.84x10"" 

Critical 

- 2.919706 

1 . 78x1 o"'' ^ , -7 . 1 3x1 0"^ 

Above critical 

- 3.019696 

1 . 63x1 0"'' ° - 6 . 78x1 0"^ 

Above critical 


Such critical values of 0, the non-dimensional 
a 

stress — ^ obtained for different values of ^ ^ are 

tabulated in Table 3,2 and represented graphically in 

a 

Fig. 3,2. It may be seen that as r 0, the value of 

Lj 

the compressive stress tends to a large value J it is 
obviously so because the motion is constrai’^ed to be in the 
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axisymmetric mode. Actually for smaller value of — , 

Lj 

the instability is likely to be in the asymmetric mode, 

TABLE 3.2 : VARIATION OF STRESS PARAMETER O = - WITH - 

aa L 


a 

L 


Pi 

Stress parameter 
' 0 

0.8 

0.94971 , 

0 

- 18.4548 

1 .0 

0.69986, 

0 

- 2.9199 

1.2 

0.52441 , 

0 

1.7586 

1.4 

0.39247, 

0 

~ 1.3642 

1.6 

0.28967, 

0 

- 1.1776 

1 .8 

0.18910, 

0 

1.0742 

2.0 

0.07437, 

0 

- 1.0111 

2.2 

0 

0.12304 

- 0.9701 

2.4 

0 

0.1 5208 

- 0.9548 



LEGEND 

<p,(r) L.H.S. OF EQ.(3-24) 
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CHAPTER ; 4 


STABILITY OF A THIN-WALLED VISCOELASTIC 
CYLINDER UNDER TORSION 


4.1 INTRODUCTION AND PROBLEM FORMULATION 

In this chapter, the critical state of stress for 
instability in a geometrically perfect viscoelastic thin 
cylindrical shell under torsion is established. Although, 
the problem of stability of geometrically perfect elastic 
and elastic-plastic thin cylindrical shells received the 
attentions of many researchers, it appears that the 
corresponding problem for a viscoelastic material has not 
yet been investigated. Here, the basic formulation developed 
by Neale (1973) for elastic-plastic cylinders has been 
extended to include the case of viscoelastic cylinders. 


A geometrically perfect thin viscoelastic cylindrical 

shell is subjected to a torque M applied at both ends. At 

any instant of time, let R , t and L be the mean radius, 

thickness and length of the shell, such that ^ < < 1 and 
R 

— < < 1 , If the boundary conditions are such that the cylindrical 

Ij 

form of the shell is preserved during its motion, the only 
nonvanishing stress components, with reference to a cylindrical 


polar coordinate (r , z) are 

M 


02 


7.Q 


2 % R'^t 


... (4.1) 


id 


all other 


0 . 
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In the current state of equilibrium, it is assumed 
that the cylindrical surfaces are free from any nominal 
traction-rate , i.e, , 

• • "t 

s = = s =0 at r = R- “- and at r = R + s— 

rr rQ rz 2 2 


Using the relation (1 «39) between s . and 
conditions (4*2) imply that 


B ^ * 

— J:J- 
t 


(4.2) 


D 


rr 


= 0 


D t 

^ ^r9 
B t 

^ rz 
B t 



a 

0z 







/ 




(4.3) 


In view of the small thickness of the cylinder, the 
preceding relations may be assumed to hold for all values 
of ; consequently 


D TT 

rr 


B 

D t 

D t 

B t 


0 , at all r 




(4,4) 


which suggests that 6^^ = 6^ = 0 , as is the case in 
thin shell analysis. 



68 


The components of the strain- rate tensor 
^ 1 / 

“ o vv, . + V. .) and those of the rotation tensor 
-) 

“ (v. _ - . ) in polar coordinate are given by 

j-j ^ » <J J f -*• 


e, = 

IG 


V dv 
r ^ dr ^ 


1 

2 oz dT 


1(1 ^ 

2 r * 60 jp 0r 
,u 1 dv V 

C“ + “ — 
r r 3 0 


1(01 + 1 Dl) 

2^5z r 80^ 


1 , 0u aw , 

2^85 ^ ^ 

1 ,8v ^ 9w, 

H — —J 

2 az r 30 


1( IJi ) 

2 5 z Q T 


1 /6 V 1 6 w 


2 9 z r d 0 


... ( 4 . 5 ) 


l(Jl lii -?2.+ 

2r30arr'^ 


1 (1^ 3 u a V ^ 

2 r 36 3 r r 


1 

2 a z 0 r 


2 ^a z r 


... ( 4 . 6 ) 


I 

I 
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In view of the stress distribution (4*1 ) and using 
the expressions (4«5) for8j_j» "the governing constitutive 
equation (1 .40) can be expressed in the following form in 
terms of the physical velocity components u^v and ws 


rr 




^tQ 


0r 


rz 


s, 


zr 

6z 

®Z0 


= (2Q+2 77D)u„ + 


S00 = (2Q + 2?7D) (“ + “Vq)+p 


— w 

r , 


e 


= (2Q + 2?7D)w ^p-Tv 


( Q + D) (1 u 0 “ “ + V ^ 


r 


1 


( Q + T) D) ( » ^ fl " - + V J + T 1 (u , - w J 


h'i 


2 fZ , 

V 


= (Q+7?D)(u +w )“'T 0 - ~ + V ) 

M 2 r I' ,r 

= ( Q * r, D) ^ T 1 (1 Z - 

= ( Q I- 7? D) (v + 1 w q) - ( “ t- “ Vo) 

,z r r r 

= ( Q )■ D) (v ^ + 1- w o) - w ^ 


,z r ",0- 


/ 


... (4.7) 


where, as "before, a comma signifies partial differentiation 

d 

with respect to the suffix that follows and D f ^ , 
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Z|,2 STABILITY CONDITION 

A general stability criterion and its derivation 
are presented in Chapter 1 , For the sake of convenient 
reference, the criterion (1,43) can be expressed as a 
requirement that the functional 


I 





] d V 


... (4.8) 


be positive for all continuous velocity fields Vj^^. 
Substitution of the rate constitutive equations (4,7) in the 
functional (4.8) leads to 




.2 


r dr d© dz 


+ /;/[ (Q + n D) 


£ 2u^ + 2( 

9 




2 

+ 2w^ 


z 


+ ( 1 u e - ; + V (V 1 w_g)2 ^ (u , + 5 


2 / 1 s 

-Tl-Wn ( -“•(•-Vq) +2v . W 

X* 2?^^ ^ ^ ^ ^ 


3 . 1 V 

— ^ 0 ^ P (v p - :” ) w } ] r dr d0 dz 

2j, f'" 2 ,r r 


• • • 


(4.9) 
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Taking into consideration the usual assumptions of the thin 
shell theory (cf,, e,g, Novozhilov 1953) the velocity 
fields are expressed as 


u = u 




V 


V + 




w 


V - ( V - U e) 





... (4.10) 


= w - a u^z J 

where u(0,2) , v (0,z) , w (&,z) are the velocity components 
at the middle surface of the cylinder along the r 0 z - 
axes respecbively, and ^ = r - R,represenbs the distance 

from the middle surface. A substitution of the velocity 
fields (4.10) into (4.9) and subsequent integration through 
the thickness of the shell leads to the following condition 
for stabilitys 


//Rtp [u + Uq + 2(1 +t)u 0 (v-U 0 ) + (•l+3t)(v-u 

9 9 , y ' 

+ ( w -H R u ^ 2t(1 + t) u „(w + R u ) 

+ (1 + 3t) u ^2 ^ 


Eq, 4.11 contd.. 
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1 + t/2R 


+ // (Q + T D) [ 2 < u + u^eg) In (' ^ ^ 


+ 2(v^q - U^gg)^ (| ) 


+ 4(u + U^gg) (v^Q - U^gg) ) 

+ 2Rt(w ^ + R u + 2R^t (1+3t) u 

IP Z I ZZ p ^ 


+ (v - U^g)2 (i ) 

+ 4Rt (v 2, ' ^^,0 ■'■ ^ ^,Gz^ ^ 

+ 2t (v^2 - ^,G ■*■ ^ ^,0z ^ 

+ Rt u^ 3 dO dz 

,z 


1 +t/2R 
1 -t/2R 


-SfT [2t{~ w^g (v^o - n^gg) ' (u +U^0 q) 

1 + *t/ 2 Rv>*^ - \ f Tt 

+ 2 In ( — — ) ( u + n gg) ( W g 

^ L /on • 


^ - t/2R 


",ez' 


+ aR^t ■" “,zz> 

+ (1+t) (V ^ ^ ^zz) 

- (1+t ) ■ ^,ZZ ^^Z ■ “,0Z^’ 

- JtUoU.Ddedz > 0 

... (4.11) 

12R^ 


in which t 
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The critical velocity distributions u,v, w are 

those which minimize the forgoing functional and, therefore, 

the EUler equations of these integral must be satisfied. 

Under torsional loading, the deformation pattern of the 

viscoelastic cylindrical shell may be assumed to consist 

of a number of circumferential wave like patterns which 

spiral around the cylinder from one end to other ; a solution 

may be sought in the following form 

W z "h 


u = A. sin ( - n 0) exp (u t) 

' R 

V = cos ( ■ — -n 0) exp (w t) \ (A. 12) 

2 R ^ 

W = A, cos ( - n 0) exp (w t) 

^ R / 


where T = 


m Tc R 


, and m, n are positive integers 


L 


5 = iw , i =^“1 , w being the frequency of motion. 

Substitution of (4.12) in (4,11) results in, 

IT L t Pq o 2 ? .2 

( ) (p WR^) [ A^ + n ilj + n A^A 2 + 

R ' 

+ (T il, + Ag)^ - 2(1 + t ) T A^ (T A^ + Aj) 

+ (1 + 3t) A^ ] 


Eq. 4,13 contd.,. 
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+ ( — ) (Q + 77 w) [ T £2(1-n^)^ A^} + £ 2n^ (nA^+A 2 )^l 

+ £4n(l-n^) ( n A^ -t A 2 ) A^ } 

+ £ 2 ^^ ('P A^ + A^)^ + 2(1 +3t) IP^A^ 

- 4(1+t) A^ A^ + A^)} 

+ £(nA^ +A2)^} 

+ £4(n A^ - + T ( n W A^ + n A^)^ 

- 2n 'P ( n A^ + A 2 ) (A^ + A^)} 

+ 'P^ A^ ] 

TC Xj *t 

- (_^ ) T [ 2 £ A^ (n A, + A^) ~ (1-n^) 

+ T £2n (1-n^) A^ (A^ + \P A^)} 

- 2n 'P^ A^ Aj + A^) 

+ 2(1 + t ) {'P^ (n Aj + A 2 ) A^+ A^)+n 'P^ A^3 

- A^ ( n A^ + A^ ) + I n A^ ] 

> 0 ... (4.13) 

in which ^ ■t/2R. 

lnr« 1 

‘l - t/2R'^ 

T = 

t 

R 
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The Euler equations of (4.13) yield three linear 
homogeneous equations in and A^, which on simplification, 

leads to 


^11 



^12 

^2 

+ 

^13 

*3 

= 0 

Or 


+ 

^22 

^2 


^23 


= 0 

^31 


+ 

^32 

^2 

+ 

^33 


= 0 



(4.14) 


where, 

+ 2n^ T 

+ 22n^ - 4n^ - 4n^ ^ ] 

- (~ ) [ 4n T (1-n^) T + 8 n t + 4n^ 'P - n ^'] 

Q 

-2 2 

P ^ R o - 9 

+ (_ ) [ 2(1 +n^) + 2t ] 

Q 

71 W T ^ _ P“^^^ 

^12 = ^ - 2n '!')-(- ) (2T^ t) + ( ) (n) 

p = (1 + ™ ) (2n^ ^ T - 4 t - 2n^T') 

1 ^ Q, 

“ ( - ) [ 2n^ + 2n ( l-n^ ) T + 2nW^ t ] 

Q 

-2 2 

p 

+ ( ) [ - 2 ( 1 + t ) ] 

Q 
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^ 2 -\ = ('1 + ~ — ) ( “ 2n “ 2n 'P + 2^^) - ( ^ ) (2^^^ ^ ) 

P S2r2 

+ ( ) (n) 

Ui 

7} W p 

^22 ~ ■'■ ) ( 4n + 10) + ( ' ■' " — ) (2) 

Ui Q 

r) 0) o o 

P = (1 + ) ( - 2 nW) - ( — ) [ an^ + 2r (1 + t )] 

^ Q 

^ p "7. p 

^31 “ — ) ( 2n 'll T - hW^ t - 2xr 

T 

- ( - ) [ 2n^ + 2n (l-n^) T + 2n t ] 

Q 

-2 2 

p 

+ ( ) ( - 2^ t ) 

P^2 = (1 + ) ( -2ii 'P ) - ( J ) [ 2n^ + 2'!'^ ( 1 + t )] 

-2 2 

r/ 03 5 _ P R 

P = (1 + — ) (4 'P‘=^ + 2n T ) + ( — ) (2) 

0 Q 

For the exisbence of a non-trivial solution of 
(4,14), the determinant of the co-efficients of the 
matrix must vanish, i,e,, 


^11 

^12 

^13 

I 


^21 

^22 

^23 

= 0 

... ( 4 . 15 ) 

^31 

^32 

^33 
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This leads to a cubic equation in non-dimensional 
stress t viscosity and frequency (after non-dimensionali** 

zation) , which has to be solved for a particular set of 
geometric and material parameters. The behaviour of the 
load-frequency relationship has to be studied by varying 
the stress VQ gradually, keeping the other parameters fixed. 
If it is found that the imaginary part of the frequency is 
positive, for a particular value of stress, the deformation 
is stable and when the imaginary part becomes negative the 
deformation is in an unstable range. The stress at which 
such a transition takes place is termed as the critical 
stress. 


4.3 NUMERICAL RESULTS AND DISCUSSIONS 


For numerical computations, the non-dimensionalised 
stress, frequency and viscosity parameters are introduced as 
follows’. 


0 » 


2Q 


P R'^ 

; Q = T ( “ — ) 
Q 






V ( PQR^ 


As the equation (4,15) is fairly involved, a trial 
and error scheme has been initially adopted to study the 
stress-requency response of the system. At a sufficiently, 
high value of stress it has been observed that the real 
part of the frequency is zero while its imaginary part is 
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negative* whereas at a lower value of stress, the real 
part of bhe frequency is non-zero and the imaginary part 
of the frequency is positive. It is therefore, natural to 
expect that at some value of stress either imaginary part 
is zero or both the real and the imaginary parts vanish. 

This value of the stress will be termed critical because 
it marksthe transition from a stable to an unstable 
configuration. It was verified at a later stage that there 
is a value of stress at which the frequency q is zero. 

Next, for a given R/t ratio, the critical non- 
dimensionalized stress parameter ® ^ ^ been obtained 

by minimizing (4,15) with respect to n and keeping the 
frequency parameter Q =0, The minimizing value of n 
is always 2 while the minimizing value of T lies betvreen 
0,2 - 1.0, higher values corresponding to larger values 
of R/t, 

Once the critical stress parameter is determined for 

a particular R/t ratio, stress is increased gradually from 

the value zero to the critical value and slightly beyond 

and the frequencies are determined for several values of 

viscosity parameter (0.01, 0,1, 0,5) keeping n and 

constant /unchanged. These are shown in Tables4.1 and 4,2 

for two values of R/t ratio (20, 50), The stress-frequency 
are 

relationships/also shown schematically in Figs. 4,1 and 4,2, 
The eleventh row in Tables 4,1 and 4.2 indicates the state of 
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the frequency parameter when the stress is just critical 
and. the twelfth row shows that the imaginary part of 
the frequency parameters are negative, when the stress is 
more than the critical value. Also, the real parts of 
the frequencies are zero at these stress levels, 

T 

The values of the critical shear stress — 

2Q 

presented in Table 4,3 and Its variation with R/t is shovm 
graphically in Fig, 4,3. The change in critical stress is 

R- D 

found to be negligible for higher values of - (L,e, — '>100). 

t t 
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TABLE 4.2 : VARIATION OF FREQUENCY WITH STRESS i ~ = 50 ; n = 2 ; ^ = 0.3752 
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TABLE A'. 3 • VARIATION OF CRITICAL STRESS WITH 

t 


1 

1 

n 

,,, mu R 

L 

0 = ( JL 

cr 2Q 

10 


0.2342 

0.1313 

20 


0.2667 

0,0516 

30 


0.2872 

0.0267 

40 


0.3124 

0,0177 

50 

2 

0.3572 

0.0128 

60 


0.4712 

0,0098 

70 


0.5455 

0,0096 

80 


0.62S1- 

0.0094 

90 


0.7713 

0.0093 

100 


0.8244 

0.0091 









t-0 = ly 


64 



ency with shear stress 




CHAPTER : 5 




CHAPTER : 5 


STABILITY OF A THIN RECTANGULAR VISCOELASTIC 
PLATE UNDER INPLANE LOADING 


5 J INTRODUCTION AND PROBLEM FORMULATION 

In this chapter, the stability of a thin rectangular 
viscoelastic plate under inplane direct loading has been 
investigated. The problem formulation is that of a plane 
stress condition. 


An incompressible homogeneous viscoelastic solid 
in a state of plane stress is considered during the process 
of continuing deformation. The current distribution of 
stress, together with other mechanical properties at this 
instant is supposed to be already determined. The current 
shape of the body is rectangular and of dimensions a x b x h. 
A set of rectangular Cartesian coordinate axes is chosen 

such that the plate, with its edges x^ = 0,a and - 0,b 
and with the plane faces =+ ” , is simply supported on 
all four edges. Whenever convenient, the coordinate x^ will 
be replaced by x,y,z and the velocity components by u,v,w, 
respectively. Referred to this frame, the instantaneous 


stress field is given by 




0 

0 


0 


0 


0 
0 
OJ 


(5.1) 
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During further infijiitesimal deformation, when the 

system is perturbed from the current state to a neighbouring 

state, it has been assumed that the ends x = 0,a and y = 0,b 

are frictionless (simply supported). Also the faces 
h 

z - + “ are assumed to be free of nominal traction rates. 
These conditions are expressed as 











0 

[ 


J 



• 1 

/ X = 0,a 

... (5,2) 




«i 3,, = 0. 0 = 2,3 J 

I 

i 





0 



J 

• 




J, y = 0, b 

... (5.3) 

T • 

D 


ss 

j “0 , 3 = 1,3 

I 


h 


- 

n. = 0 

h 

- ± 2 

... (5.4) 

where 

^i 

is 

the unit outward normal to the boundary 


surfaces. 






For 

the present case, the 

rate const! butive 

equation 


for the viscoelastic material, as described by (1.40), can 
be expressed in the following explicit form 


'xx 

= (2Q -i- 2r) 0 - a^) u^ 

+ P 

yy 

= (2Q + 2t) D - a^) 

+ P 

zz 

= (20 + 2 rj D) w 

+ p 


Eq, 5.5 contd.. 
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s 


pc 


(Q + T) D - 
(Q + 77 D- 



+ V ) 



V. 



+ v ) + Qf 


2 


u 


y 


# 


S 


xz 


s 


zx 


s 


yz 

zy 


(Q + D - ^ ) ( 

0 

(Q + rj D - ) (u + w„) 

2 z X 

(Q + r) D - ^ ) ( + Wy) + Wy 

(Q + Tj D - ) (v + w ) 

2 2 y 

... (5.5) 


where, D = , and a subscript denotes partial 

9 1 

differentiation. Also, the velocity is replaced by its 
physical components u,v,w. It is to be noted that these 
velocity components are not independent but are connected 
through the incompressibility condition 


u + V + w 
X y z 


0 


... (5.6) 


The boundary conditions (5,4) at z = + — suggest that 


zz 


(2Q + 2T7D) w + p 


« 0 




zx 


= (Q + 77D-2i)(u + w ) = 0 


^zy " (Q + ^ D - — ) ( Vj^ + Wy) = 0 




\ ... (5.7) 


The first equation of (5.7) can be satisfied by imposing 

that s is actually zero through the plate thickness, 
z z 
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The remaining two conditions in (5,7) are automatically 
satisfied by adopting the usual thin-plate assumptions that 
6 and 6^2 much smaller (and can be taken as zero) 
compared to in-plane strain components. Hence, the constitutive 
equations (5,5) take the following simplified form 


®xx = ( 2Q + D - CT^) + p 

Syy = ( 2Q + 2q D - + p 

S 22 = ( 2Q + 2h D) w^ + p 

V = ( Q + n D - H. cr^ 

*^2 

Syx = ( Q + ^ D ( "^y 

Sxz = '^1 

s =0 
zx 

®yz " ^^2 ^y 

s = 0 

zy 


...(5.8) 
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5.2 STABILITY CONDITION 

As in the previous chapter, the criterion (1 ,43) 
is used for investigating the stability of the system. 
Substitution of the constitutive equations (5,8) , along 
with the incompressibility condition (5,6) in (1.45) leads 

to the following condition for stability! 

, h 

P I I ft I- ^ V + ] dx dy dz 

0 0 

a b - 

^ f S fl [ (20 4 

0 0 

2 

+ ( 2 Q + 2 r) D - + ( 2Q + 2r) D) 

+ ( 0 + r,0 - 2) ^ 

+ ( 0^ vj + °2 ) + ( '^^ + °2 ^ 

> 0 •••(5*9) 

An approximate velocity field is chosen which 
incorporates the usual engineering assumption that line 
elements perpendicular to the middle surface incipienbly 
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remain straiglii: and normal to the middle surface j its 
familier form is 


u = - 2 w. 




X 


V « - z w. 


W = W (x,5r,t) 

Substituting the velocity field (5,1 O) 
and subsequent integration through the 


p n 

12 

a 

/ 

0 

b 

/ 

0 

\ — 1 

.2 

+ w + 

y 

12 

W J 

2Q 

h3 

a 

b 


77 D 

a 



/ 

S [ 

( 1 + 

- ^ 

12 


0 

0 


2Q 




277 D 

a 

+ 



+ 

(2 + 


- J_, 

2 





Q 

2Q 

' xy 


a 

b 



a 


+ 2Q h 

f 

0 

/ 

0 

2Q 

2 

w^ + 

2 

2Q 

I 1 

> 


... (5.10) 

criterion 
in the stability /( 5. 9), 

thickness, results in 
dx dy 




n D 
Q 


a 

”” )w^ 
2Q y-y 


] dx dy 


dx dy > 0 


... (5.11) 

Let the solution of w be taken in the following form 


w = sin vx sin ky exp ( 5 t) ,,, (5.12) 

m 7t 

where, = — ~ k = , m = i m , m, n being 

a b 

positive integers, co is the frequency and i It is 

expected that this field will furnish a 'close* upper bound 
for the critical stress. Use of (5,12) in (5,11 ) results 
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-^2 


in the following requirement for stability 

p 5 
2Q 
12 


r' OO 0-1 

[ ) + ( ^ + ^ ] 




Q 


[ { 


- ( + K^)) ( ^ k^) ] 


h'= 


2Q 


2Q 


> 0 

which on simplification leads to 

a 


2Q 


p/- 


2 

2Q 

Tr- 


1 - s' 


1- S 
where, 


p'. r?(^f 1- C1 . 


12 


~ 


r'- 'fr (1 + -r- ^ 


2 ^ ^ 


<1 


... (5.13) 


... (5.14) 


^ , mb na ^2 P (jJ a 1 c t \ 

) ( — ) + (1+S ) 

na mb 2Q 


mb na 2 P m^a^ 


a' 12 


Q 


na mb 


2Q 


.0 


n % 


.(1+S'> 


11^ r 2/ 2 ?,h 0 -\ 

and S’ « [ m (—* ) + n (t^) ] 

12 a ® 


The inequality (5.13), after substituting (o = i m, 
can be rewritten in the following form 

^2 _^3(ij^C<0 ••* (5,15) 
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in which 

A P / 2 2 ^2 

A=_(v + k- + ~ ) ... (5.15a) 

B = ~~ + k^)2 (5.15b) 

C = [ I ~ “(v^+k^)} (^ v'^ + ^k^) + (v2+k^)2] 

... (5.15c) 

It is to be noted thet A,B,C are real quantities while w in 
general, is complex. 

The inequality (5.14) further leads to 

(o) - ) ( 0) - oj^) <0 ... (5.16) 

where, 

B 1 P 1/2 

^1,2==^~ i _ [ ~ B + 4AC] ... (5.17) 

Evidently, the first term in (5.17) represents the imaginary 
part of the frequency (and is normally positive when the 
motion is stable), while the second term in (5,17) represents 
the real part of the frequency if the expression within the 
square bracket is positive* Depending on the values of A,B,C, 
it is possible that the expression ( *- + 4AC) in (5.17) 

may be negative. In that case, the real part of the frequency 
is zero and the imaginary part of the frequency is non-zero. 
Alternatively, then, 
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Hence, for imaginary part of bhe frequency to be negative 

2A 2A ^ 

which yields 4AC < o. 

Since, A > 0, for imaginary part of the frequency to * be 
negative, 


^ ^ ^ (5.19) 

Using (5,15c), the preceding condition implies that 


[ 


12 

i7 


- ( 


n^ 71^ 


ir ^ <7 2 2 

^ n IX 


Qr 


^ TT 


■ 2Q a2 ‘ 2Q ^.2 


] 


r n2 ,2 

] < 0 


which after simplification, can be rewritten as, 
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2Q 


2 /h ^2 
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2Q 


ra^(— )^ E* 
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whe re , 


•n;' 
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mb 


E = - 1-2 , . ...na 


na ■ 

— V 
jnkJ. 


[ 1 - ( - )' 
12 a 




. .. (5.20) 


Hence, for imaginary part of the frequency to be negative 
(i.e, when the motion is unstable), the condition (5,20) must 
be satisfied. Further, as evident from (5,18) and the 
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foregoing explanations whenever imaginary part of the 
frequency is negabive, its real part is zero. Obviously, 
when the motion of the system passes from a stable state to ^ 
an unstable one, due to Increase in the stress ( and/or ), 

2Q 

the frequency passes through zero. Hence the critical value 
of stress would correspond to the situation when m f 0. 

It is also unlikely that the imaginary part of the frequency 
will attain a negative value some where inbetween during 
gradual increase oC stresses from zero to the value obtained 
From (5.20) . 


5,3 numerical results and discussions 


For numerical computations, the stresses, frequency 
and viscosity of the material are non-dimensionalised as 
follows 5 


0 = -L ; 0 

' 2Q ^ 


Q =i 


pa 


2Q 


= V( — ) . « ; = 




... (5,21) 


The variables used in the numerical computations are 
The plate dimensions characterized by — (i.e, ratio 
of the thickness of the plate to its side in the 
y - direction) and aspect ratio - . The value of 

is taken as 0.01 , while the ratio | is varied. 

b 
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ii) Material property is characterized by non-dimensional 

viscosity parameter (as defined in Bq. 5.21)* 

Three values of X^ ( =0.01 , 0*10 and 0.50) are 

baken to sbudy its effect on the frequency parameter 0* 

ili) 'Dwo particular load-conditions are considered 

a) ^2 *= 0 plate is loaded in the 

X“direction only 

b) 9^ / 0 , 02 0 • ibe plate is loaded biaxlally. 


When tho plate is loaded in one direction, the 


critical value of bho stress 0^, as obtained from (5.20) is 

, Jii f fi H. ^ ,2 

. 2 /h ^2 12' na mb 


0 


1 


= - ) 


na 


[1 -"_ Cm 2 (^)^ h .„2 ( 1 ) 2 !] 


(5.22) 


a 


Tho above expression is minimum for n = 1 and various 
values of stress parameter 0^ are plotted against the aspect 
ratio ~ for different values of m (v/ith n=1 ) in Fig. 5,1 ♦ 

Tho lower envelop of bhese curves, as indicated by the 
thicker lino, denotes the minimum value of the stress parameter 
0-j below which the plate is certainly stable, 

Whon the plate is loaded biaxially critical value 
of tho stresses 0^ and ©g sire determined from (5,14) which 
represents a system of straight lines in 0^ and ©2 with 

"n / D ^ 

intercepts of r and ^be abscissa and the 

) (1 - s’) 

ordinate respectively, Tho system of straight lines 
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in — o and — ^ are first plotted with different 

(| )2 ( ^)2 

m and n (for a particular aspect ratio) and the critical 
envelop is determined. For various aspect ratio , the 
cribical envelops are shown in Figs. 5.2 and ’5,3 « The following 
observations can he made from these two figures: 

b 

1 , Once the ratio ~ (or ^ ) and — are fixed, 

it is easy to determine the stress parameters ©2 
and 0>j which will cause instability. 

2 , It is obvious from Figs, 5,2 and 5,5 that the 

instability is impossible when parameters 0^ and 
aro both tensile. But the possibility does remain 
when one of the stress parameters is tensile while 
bho other one is compressive, as well as when both 
stress parameters are compressive. 

Application of tensile stress (load) in one 
dirocbion, increases the value of the critical 
3 Irons, 

Onco the cribical stress is determined for a particular 
value of - , ~ and ^ > "the stress parameter (and hence the 

load) is iLroaLd grad^llly from the value zero to its critical 
value and slightly beyond (say, in steps of ^ The frequency 
parameter o is determined for three values of viscosity 
paramoter (0.01 , 0.10 and 0.50). These are tabulated in 
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Tab'Jes5,1 bo 5,3« The variation of real and imaginary 
parts of the frequency parameter with stress parameter, is 
shown schematically in Figs, 5,4 to 5,6, The eleventh row 
in Table 5,1 - 5,3 indicates the state of frequency parameter 
when the stress is just critical and the twelfth row shows 
that the imaginary parts of the frequency-parameters are 
negative when stress is more than the critical value. Also 
the real parts of the frequency-parameters are zero at these 
stress levels, 
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CHAPTER : 6 

CONCLUSIONS 


The pfesent investigation has been concerned with 
some stability problems in a finitely deformed incompressible 
linear viscoelastic solid. Conditions for instability have 
been obtained in all the problems studied, 

Exhaustive formulation has been presented to 
Invo stiffs bo instability in a rectangular solid in plane 
strain. In this case, the numerical solution of the 
instability condition shows that the only possible mode of 
instability is the antisymmetric one for all possible 
material dimensions. The analytical expressions derived 
for instability in the symmetric mode yield higher value 
of the critical stress than those obtained from the 
corresponding expressions in the antisymmetric mode and 
arc, therefore, of only academic interest. The results 
obtained for the limiting case of the surface instability 
in a viscous medium have been found identical to those of 
Biot (1965). 

The condition for axisymmetric instability has been 
derived for an incompressible viscoelastic solid cylinder 
of f ini bo dimensions under axial loading. As resulting 
condition is quite involved, it has not been possible to 
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rule out instability under tensile loading. However^ the 
numerical solution of the condition yields only the 
compressive values of critical stress for different cylinder 
dimensions, 

For a long thin viscoelastic cylindrical shell the 
value of the critical shear stress (and hence the torque) 
has been obtained at which the shell becomes unstable. These 
havo boon plotted for different radius-thickness ratios, 

Tn the caso of thin rectangular viscoelastic plates, 
the conditions for instability have been examined both under 
axial and biaxial loading. The instability due to unidirectional 
in-plano tensile loading has been ruled out. Critical stresses 
have been obtained for the axial and the biaxial loadings for 
various aspect ratios. In case of biaxial loading, the 
interaction curves have been plotted to illustrate the region 
of stabilj by/ instability. 

It has been observed in all the problems that, as the 
sbross/load gradually approaches its critical value and goes 
slightly beyond, the frequency of the system passes through 
tho value zero. In other words, when the stress/load is 
below the critical value, the real part of the frequency is 
non-zero while its imaginary part is positive. The stage at 
which both the real and imaginary parts of the frequency are 
zero is termed as tho critical one because it marks the 
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bransition from a stable state to an unstable state* When 
bhe sbrcss/load attains a value higher than the critical 
one, the real parb of the frequency remains zero while its 
imaginary part becomes negative. This indicates divergent 
typo of motion and instability of the system. 

It is obvious that the loss of stability is of 
static type. This gives an impression that the viscosity 
parameter has no effect. But, this is not so. The critical 
stress Is In terms of a material parameter Q which actually 
is n function of the current stress. In fact, the critical 
stress must be expressed in terms of the material parameter 
Qq in tVic; underforiftod configuration, The relationship 
botwoon Q and depends upon the manner in which the 
body has boon finitely deformed (see, e.g, Biot, 1965) and 
It, naturally, is dependent on viscosity. 
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